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MIXING CHARACTERISTICS DOWNSTREAM OF CORE REGION OF HIGH- 
TEMPERATURE AXISYMMETRIC JETS EXHAUSTING INTO 
TRANSONIC AND SUPERSONIC STREAMS* 

By Ronald G. Huff and Kaleel L. Abdalla 


SUMMARY 

A series of rocket motors with varying exit to throat area ratios 
was tested in the 8- hy 6-foot wind tunnel to determine the effects of 
mixing on jet diameter and temperature decay at large distances (x/d > 30) 
from the nozzle exit. An approximate method to account for effects of the 
initial expansion was evolved. It was determined that the combustion 
efficiency has an important effect on jet spreading, since the unburned 
products can bum downstream of the nozzle. The data showed considerable 
scatter j however, mixing rates were, in general, lower than those ob- 
served for subsonic jets. Data for angles of attack of 5° and 10° are 
also presented, giving the respective centerline shift and temperature 
decay as a function of axial distance. 


INTRODUCTION 

Supersonic jet-spreading rates and temperature distributions at 
large distances downstream of the nozzle exit are of importance to air- 
craft an d missile designers for problems such an jet infrared radiation 
nnri engine inlet ingestion of exhaust products from airborne rockets. 
"Large distances downstream" means distances on the order of those re- 
quired for elimination of the central potential core in the jet. The 
shear region extends across the entire jet, and mixing profiles at suc- 
cessive stations are essentially similar. 

The general case is complicated by the fact that the jet issues 
initially from a nozzle having a nonzero exit angle into a region of dif- 
ferent static pressure. Studies of Jet spreading in the region near the 
nozzle exit (see refs. 1 and 2, e.g. ) have shown that the process of 
pressure- and flow-direction equalization can have a much greater effect 
on initial jet-spreading rates than the mixing process. 
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The size of the jet and the spreading rate at large downstream dis- 
tances will then depend on two independent processes: first, the 

pressure -direction equalization that rapidly changes jet size near the 
nozzle, and second, the more gradual mixing of the equalized jet and the 
surrounding stream at constant pressure. 

Available data on spreading rates near the nozzle show that the 
problem is complicated, and the data are not sufficient to obtain effec- 
tive jet size after equalization as a function of body-nozzle geometry, 
jet pressure ratio, stream Mach number, and so forth. Suitable informa- 
tion on mixing is likewise lacking. Although considerable theoretical 
and experimental work has been done with subsonic jets and streams (see, 
e.g. , refs. 3 and 4), supersonic spreading rates as a function of veloc- 
ity and density ratios between jet and stream are not available. 

The primary objectives of the present investigation were to obtain 
spreading and temperature distribution data for hot supersonic jets ex- 
hausting into high-speed streams (Mach 0.8 to 2.0) and to show the effect 
of the initial expansion on these downstream spreading characteristics. 

A secondary objective was to develop a simple method that, with the use 
of existing information on free jet expansion, would provide approximate 
answers for spreading rates. 

Three small rocket motors having a nominal thrust of 250 pounds and 
of varying exit to throat area ratio were used to produce the hot super- 
sonic jets. The motors used JP-LOX at a stoichiometric fuel-oxygen ratio 
and were tested in the 8- by 6-foot supersonic wind tunnel at stream Mach 
numbers from 0.8 to 2.0. The basic data were obtained at an angle of 
attack of zero) angles of 5° and 10° were investigated briefly. A total- 
temperature rake was used to survey the flow. 


SYMBOLS 

A area 

c* characteristic velocity 

d diameter 

dQ g "half” diameter of the jet) diameter at which temperature 

equals (T m + Tq)/2 

F,G Mach number functions 

g acceleration due to gravity, 32.174 ft/sec^ 

M Mach number 
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^max 
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r 

e 

A 

^adj 

^eff 
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T adj 


molecular weight 
oxidant-fuel mass ratio 
total pressure 
static pressure 
universal gas constant 
Reynolds number per foot 
radius 

total temperature, °R 

maximum total temperature in jet at a given station 
velocity 

total flow rate into rocket motor (LOX, JP-4 weight flows) 
axial distance from nozzle exit 

local jet centerline displacement due to angle of attack 
angle of attack, deg 
ratio of specific heats 
combustion-efficiency parameter, T c /T^ 


(Tg - T C )/(T’/T C ) = (t - l)/[l/(T) - 1)] 

ratio of free-stream velocity to equalized jet velocity, Vq/V 2 
value of A e pf for which temperature and spreading data agree 
effective low-speed value from fig- 5 
density 


afterburning parameter, 


Tc 


Tc 


Ti T c 

value of t for which temperature and spreading data agree 
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Subscripts: 

^ jet centerline, point in jet where maximum temperatures exist 

c combustion chamber 

e nozzle -exit station 

max maximum 

min minimum 

t rocket throat station 

0 free-stream conditions 

1 station at which jet is fully e:q?anded 

2 station at which burning ceases and mixing begins 
Superscript: 

' denotes theoretical combustion-chamber conditions 


APPARATUS AND PROCEDURE 

The rocket motors used were water-cooled and utilized JP-4 and LOX 
as the rocket propellant. A sketch of motor and model details is shown 
in figure 1. The three nozzle -afterbody combinations that were tested 
are sketched in figure 2. All nozzles were conical with a 15° half- 
angle and had a throat diameter of 0.5 inch. The area ratios, A e /A^., 

were 3. 44, 6. 25, and 18. 31. The base size and the boattail angle were 
the minimum values consistent with the motor structural and cooling re- 
quirements. The motors were mounted in tie supersonic test section of 
the NASA Lewis 8- by 6-foot transonic -supersonic wind tunnel. Angles of 
attack up to 10° were possible and could De obtained without changing 
the location of the nozzle exit. 

Figure 3 shows a sketch of the total -temperature survey rake and 
its mounting mechanism as installed in the transonic portion of the tun- 
nel. The sting support strut is shown in the fore position. It was used 
as a steady rest and was moved to a position aft of the sting mounting 
strut for far-downstream rake locations. For each test the rake was 
traversed as far upstream as consistent with the maximum permissible rake 
temperature of 2000° F. The rake was so alined as to measure the upper 
half of the jet profile at zero angle of attack. Rake thermocouples were 
spaced at 0.25 inch over the lower 4 inches (near the jet centerline); 
for the remainder, spacing was 1. 0 inch. 
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Propellant flow rates were adjusted to obtain chamber pressures 
(nominal) of 300, 400, and 480 pounds per square inch. L0X to JP-4 ratio 
was held at stoichiometric (0/F = 3.4). Nominal tunnel test conditions 
at which each motor was run are given in the following table: 


Mach 

number. 

Mo 

Static press. , 

P 0' 

lb/sq in. ahs 

Total press. , 

p o> 

lb/sq in. abs 

Total 
temp. , 

V 

°R 

Reynolds 
number 
per ft, 
Re 0 

0.8 

9. 48 

14. 45 

586 

3. 6 4X10 6 

1.5 

5.11 

18. 75 

633 



4. 65X10® 

2.0 

3. 15 

24. 65 

674 

4. 78X10 6 


Instrumentation and Data-Re cording Procedure 

Rocket chamber pressure P c was recorded on a pressure transducer. 
The fuel and L0X weight flows were recorded using turbine-type flowmeters. 
Chromel-Alumel thermocouples were used on the jet temperature survey rake 
a.nd were of the shielded aspirating type (see fig. 3). One thermocouple 
of this type was placed on each side of the rake at approximately the 
horizontal jet centerline. These showed fairly good rake alinement. 

Chamber pressure, JP-4 and L0X weight flows, and one scan of the 
rake temperatures were recorded as soon as the rocket motor had reached 
steady-state operation. Then, as the rake was moved upstream in steps, 
only the rake temperatures were recorded. 


Data Reduction 


Rocket motor parameters were obtained as follows. The theoretical 
characteristic velocity was computed from 



where y' T^, and Jt' c were obtained from reference 5 for the proper 
values of P c and 0/F. The actual characteristic velocity was computed 
from 
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Combustion efficiency and chamber temperature could then be obtained 

from 

r\ = [c*/(c*)'] £ (3) 

and 

T c = r\T' c (4) 

The values of P c /Pe a-nd- Te were obtained by plotting the tab- 
ulated values of reference 5 as a function of area ratio A e /A^. The 
value of the exit Mach number was then obtained from 



The measured rake voltages were converted to temperatures, and these 
temperatures were used directly to compute the temperature ratio 
(T - TQ)/(T max “ Tq). (Calibrations of similar thermocouples have shown 

recovery factors of 0, 99 -k ) 

The half diameter 5 was calculated as twice the distance be- 
tween the thermocouple having the maximum temperature and that having a 
temperature equal to (T max + Tq)/2. 

STATEMENT OF PROBLEM AND METHOD OF ATTACK 
Equalization 

Although the object of the investigation is to study fully developed 
supersonic mixing rates, the pressure- and flow-direction process that 
occurs near the nozzle exit is of direct concern, since it can result in 
large changes in the size of the jet and hence in the scale of the mixing 
process. The complexities of equalization are illustrated in figure 4(a). 

The jet and the stream first expand into base region in which the 
pressure is usually below ambient. Both streams next pass through a 
series of shocks and expansions that' simultaneously turn the flow toward 
the axial direction and adjust the static pressure to ambient. Mixing 
starts in the base region. The initial boundary-layer velocity profiles 
are altered first by the expansion at the base, next by the mixing with 
the semidead air in the base region, and again as the flow passes through 
the trailing shock and the succeeding expansions and shocks. An added 
factor is the possibility of burning within the jet downstream of the 
nozzle exit. 
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In the region close to the nozzle exit, jet-spreading rates are then 
a complicated function of the equalization process and the mixing along 
the boundary. Fortunately, if the shocks and expansions are strong, the 
process is essentially complete within 2 to 3 diameters of the exit (see, 
e.g. , the schlieren photographs of ref. l), so that the major part* of the 
mixing process must occur essentially at constant pressure. Therefore, 
if an estimate of the effective diameter of the jet at the completion of 
pressure- and flow-direction equalization can he obtained, it may be 
possible to correlate spreading data obtained at different jet pressure 
ratios and/or with different nozzles and, in addition, to correlate the 
present data with existing low-speed data of reference 4. 


Idealized Model 

The actual equalization process is far too complicated to generalize $ 
however, for the idealized case of figure 4(b) the effective diameter at 
the start of mixing can be obtained. 

For simplicity, the equalization process is taken in two steps. In 
the first, the pressure-direction equalization region (between stations 
e and l), the jet expands freely to the ambient pressure p q. Inter- 

actions between jet and stream as well as mixing and heat transfer along 
the boundary are neglected. The second region (between stations 1 and 2) 
is included in the idealized process so that any burning that occurs 
downstream of the nozzle exit can be accounted for. The afterburning is 
assumed to occur at ambient pressure. Downstream of station 2 the two 
■uniform parallel streams start to mix at constant pressure. The region 
of interest in the present investigation is, of course, downstream of 
station 2 where the mixing has progressed to the jet center and the 
constant-velocity core has been eliminated. 

The procedure used to calculate quantities at stations 1 and 2 was 
first to obtain the diameter dq from the characteristic solution of 
reference 6. (Details are in the appendixj results for each run are In 
table I. ) The Mach number Mq and the stream-jet velocity ratio A 

were then calculated using jet mass -flow continuity. 

The proper value of the afterburning temperature ratio was, of 
course, not known beforehand. However, for each assumed afterburning 
ratio, T = T£/T c , the diameter dg could be calculated using the one- 
dimensional continuity and momentum relations. (As shown in the appendix, 
the velocity ratio is unchanged by afterburning. ) Several values were 
assumed, and the proper value was obtained from data analysis, as will 
be discussed more fully after the data are presented. 
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Low-Speed Data 

The results of reference 4 on the constant -pres sure mixing of coaxial 
streams include velocity ratios in the range of interest for the present 
investigation and were therefore used for comparison. Actually the tests 
were run with cold jets (Tq/T c « l), and the spreading characteristics 

were obtained by measuring the rate at which helium in the jet mixed with 
the stream. However, as discussed in reference 4, experimental evidence 
indicates that mass and temperature spread at essentially equal rates for 
turbulent mixing. 

Curves obtained by cross -plotting the data of reference 4 and assum- 
ing that temperature ratio could replace helium concentration are pre- 
sented in figure 5. Also included is a sketch illustrating the termi- 
nology. Jet parameters at the start of mixing are d 2 , T 2 , and V 2 . At 
any axial station, is the maximum local temperature and d 0 . 5 is 

defined as the diameter for which the temperature is above ambient (or 
below T raax ) by half of the local maximum increment. Since the best fit 

to the temperature profile data is a cosine curve, it can be seen that 
do. 5 is also tk e "half diameter” of the jet. 

In the section that follows, the present data are compared with the 
curves of figure 5. The effects of stream-jet density ratio and the ab- 
solute magnitude of velocity (compressibility effects) should show up 
as discrepancies. 


RESULTS AND DISCUSSION 
Zero Angle of Attack 

Temperature profiles . - The temperature profiles were normalized by 
plotting (T - T 0 )/(T max - Tq) against r/r 0# 5 , where r q § is the "half” 

radius. The shaded area representing data for all nozzles and over the 
pressure-ratio and axial-distance ranges tested is compared in figure 6 
with an empirical cosine relation that was found to give good agreement 
in reference 4. It can be seen that the cosine variation falls close to 
the center of the shaded region except near the outer edge of the mixing 
region. 

Jet spreading. - For the spreading data of figure 7, both the half 
diameter dQ 5 and the axial distance x are expressed as ratios of the 

calculated diameter after equalization d 2 . As has been explained, the 
use of d 2 is an attempt to account for the initial expansion. Zero 
afterburning was assumed (t = l), so that d 2 = d-^. The pertinent 
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parameters for each run are listed. (Additional information is listed 
in table I. ) Also included are the curves of figure 5. 

In figure 7 , many of the data sets have slopes that are in good 
agreement with the curves for the low-speed data. A large number have 
near zero slope , however, and no reasonable explanation for this result 
is apparent. In general, if the present data are compared with the curve 
having the same velocity ratio, the present data indicate lower spreading 
rates. However, if, as has been mentioned previously, density differences 
between jet and stream affect mixing rates, the present data and the 
curves should not be compared at the same A. 

Temperature decay . - An additional possibility for the discrepancies 
becomes evident from the maximum (centerline) temperature distributions 
of figure 8. The temperature ratio (W - T o)/( T 2 " T 0) is Plowed as 

a function of x/d^. assuming no afterburning (Tg = T c and dg = d-^). 

Comparison with the low-speed curves again indicates either lower average 
spreading rates (i.e., higher temperatures) or the need for different 
A values. An Important additional point, however*, is the fact that the 
temperature ratio exceeds unity by large amounts for the low-combustion - 
efficiency runs. It is clear that the unburned fuel must burn downstream 
of the throat and that as a result of the afterburning the effective tem- 
perature at the start of mixing T2 must In some cases be greater than 
the combustion temperature T c . 

Afterburning and velocity ratios . - Three possibilities for the dis- 
crepancies between the low-speed, constant-density data and those of the 
present investigation have been mentioned: 

a) Afterburning, which changes the effective temperature and jet 
diameter at the start of mixing 

(2) Density differences between jet and stream, which can change 
mixing rates 

(3) Differences in velocity magnitudes 

An additional possibility, of course, is that the computation for 
the effective diameter after expansion is not sufficiently sophisticated. 
However, for runs 9 and 10 the jet exit pressure ratio P e /Po c l° se 

to unity (see table I), and It is probable that the diameter ratio h-^/d e 

cannot be greatly in error. For runs 9 and 10, then, one or more of the 
three preceding faults must exist. 

The data of figure 8 show that afterburning is more than a possi- 
bility; the problem is that the proper amount of afterburning is not 
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known a priori. The problem was handled by using a technique in which 
an afterburning ratio (t = T 2 /T c ) was assumed. The diameter d 2 was 

then computed; the quantities x/dg, (3 Lq # g/dg, and (^max “ T oV( T 2 “ T 0^ 

were obtained; and the data were replotted. The results were superim- 
posed on the curves of figure 5, and effective low-speed velocity ratios 
were read, off for both the spreading data and the temperature data. The 
process was then repeated for a different assumed afterburning ratio. 

Typical results of the procedure are presented in figure 9. Effec- 
tive low-speed velocity ratio 7\ e ff is plotted as a function of assumed 

afterburning ratio for the data of run 10. It can be seen that, as after- 
burning ratio increases, the effective velocity ratio increases for the 
spreading data and decreases for the temperature data. Moreover, the 
afterburning ratio for which A e ff equals the actual A is different 

for the temperature data (t = 2.9) and for the spreading data (t < 1. 0). 
There does exist, however, an afterburning ratio for which both data 
sets are satisfied. For run 10 the "adjusted" values are T adj = 2* 46 
and Ag^j = 0.44, and the velocity-shift ratio Aq^j/A is 1.91. 

A summary of adjusted values is included in table I and is also 
shown in figure 10, where both ratios are plotted as functions of com- 
bustion efficiency. The symbols indicate stream Mach number (symbol 
shape), jet exit pressure ratio (tail location), and nozzle area ratio 
(solid or open). Also included are curves of constant 

9 ~ (Tg - T C )/(T' - T c ). The numerator is the afterburning temperature 

increase; for stoichiometric o/F (3.403), the denominator is a measure 
of the maximum possible increase. The value of 9 is therefore essen- 
tially the ratio of actual to maximum temperature rise for all data ex- 
cept that of run 7 (o/F = 2.73). 

The velocity-shift-ratio data show a scatter that has no definite 
trend with either stream Mach number, pressure ratio, area ratio, or 
combustion efficiency. Most points fall within Aq^j/A = 1. 5±0. 3, thus 
indicating definitely lower mixing rates than would have been determined 
using the calculated velocity ratio A and figure 5. The afterburning 
results show that, independent of nozzle type and operating conditions, 
the temperature rise is about half of the maximum value. 

Using an uncertainty in velocity-shift ratio of ±0.3 with the curves 
of figure 5, the situation in supersonic hot jet mixing can be summarized 
as follows: With the simple approximation for the effect of initial ex- 

pansion and using a velocity-shift ratio of 1.5 on the curves of figure 
5, the resulting jet size and temperature will be within ±18 percent for 
x/d^ less than 100. Whether density or compressibility effects are 

responsible for the lower mixing rates is not known. 
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Angle of Attack 

Because of the limited amount of angle -of -attack data obtained, ex- 
tensive analysis was not made. The discussion is limited to presentation 
of the results and observations of the major trends. 

Centerline shift. - Shift variations with angle of attack, jet pres- 
sure ratio, and stream Mach number are plotted in figure 11. Centerline 
shift is expressed as a fraction of the effective diameter dp From 

figure 11(a), it can be seen that increasing Mach number decreases the 
offset ratio despite the fact that a pres sure -ratio increase accompanies 
the Mach number increase. For the larger-area-ratio nozzle, however, 
offset ratio is almost independent of Mach number (compare figs. ll(b), 
(c), and (d)). For the range of variables investigated, shift ratios 
were approximately 1 for 5° angle of attack and 2 for 10°. 

Mixing properties. - In figures 12 and 13 are shown the effects of 
angle of attack on mixing characteristics. It is apparent from figure 12 
that the spread of the jet, in the vertical test plane (no measurements 
of jet spreading were taken in the plane normal to the vertical plane), 
is not appreciably affected by angle of attack. On the other hand, fig- 
ure 13 shows that angle of attack definitely decreases the jet maximum 
temperature. 

Unfortunately, most of the data in figures 12 and 13 are for over- 
expanded nozzles (p e /po < 1. 0). However, the data in figures 12(a) and 

13(a) at Mach numbers of 1.5 and 2.0 are for underexpanded nozzles. These 
data agree with the preceding conclusion. It should be noted, however, 
that the differences of jet maximum temperature with angle of attack are 
smaller if the nozzle is overexpanded. 


SUMMARY OF RESULTS 

The following conclusions on downstream mixing are drawn from the 
results of a test conducted in the NASA Lewis 8- by 6 -foot wind tunnel 
using JP-4 and LOX as propellants in small rocket motors operated at 
various combustion efficiencies: 

1. The jet temperature profiles at various axial distances were 
similar and approximated those obtained in existing subsonic data. 

2. Mixing rates were lower than those for the low-speed cold jet 
data of reference 4. 

3. Using a simplified method to account for the initial expansion, 
for afterburning, and for the lower mixing rate, the spreading and tem- 
perature decay data showed a scatter of ±18 percent when compared with 
existing low-speed data. 
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4, Because of afterburning and the resulting change in effective 
jet size and temperature, combustion efficiency is an important parameter. 

5* Angles of attack of 5° and 10° shift the jet centerline by approxi- 
mately 1 and 2 jet diameters, respectively, at Mach 1.5. 'While the jet 
maximum temperature decreases with angle of attack, the spread of the jet 
remains nearly constant. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 13, 19*9 
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APPENDIX - METHOD USED TO CALCULATE JET CONDITIONS 

AT STATION WHERE JET AFTERBURNING CEASES 

Mach number functions that will be useful in the derivations that 
follow are: 

1 

M^l + M 2 

and 

o-Vr -; -v 

yi + I ~ 1 m 2 

Charts of F and G as functions of M were prepared for several 
values of y. 

From flow continuity between the nozzle throat and station 1 and 
assuming T^ = T c and p^_ = p q, 






In the area ratio A^/A^ A^ was assumed to be the maximum area in a 
free expansion. Numerical values were obtained as functions of P c /pq 

and M e from the characteristic solutions of reference 6. The bracketed 
term on the right in equation (Al) was neglected, since, for the maximum 
pressure ratio tested (P c /pq = 150), it differs from unity by less than 

1 percent. Mach number was found from equation (Al) and the curves 

of F against M. 
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For the assumed const ant -pres sure afterburning between stations 1 
and 2 (pq = P£ = Pq)j equations can be written for conditions at station 
2 as follows: 

Continuity: PqVqAq = ppV-Ag 
Momentum: (piVqAq)Vq = ( p 2 ^ 2 ^ 2^2 
Obviously, 


Vi = v 2 


From equation (A2) and the definition of G, 


(A2) 


^2 



(A3) 


The Mach number M 2 can then be found from the G against M curves 
for any assumed value of afterburning temperature ratio t = T 2 /T c . 


Finally, the jet diameter at station 2 can be found by rewriting 
the momentum equation: 

n M f A i = * 2 M I^2 


h % E 

di - M 2 y r2 


(A4) 


The ratio of the stream to jet velocity 7\ at stations 1 or 2 was 
obtained from 


_ _0 _ L % 

= v x = y Tc )u 0 Gl 


(A5) 
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ABLE I. - ROCKET -MOTOR OPERATING CONDITIONS AND PARAMETERS 
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Figure 1. - Typical rocket motor. 





Aspirating thermocouple 



.nstallation in 8- ty 6-foot wind tunnel. 




(b) Idealized equalization. 


Figure 4. - Jet-stream pressure- and flov-direction equalization. 
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Figure 5. - Low- speed spreading curves (from crosB plots of concentration 
data of ref. 4). 







Figure 6 . - Temperature profile . 
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Figure 8. - Temperature decay. No afterburning correction 
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(c) Area ratio, A 0 /A t , 18.31; stream Mach number, 1.5. 
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(d) Area ratio, Ag/Afc, 18.31; stream Mach number, 2.0. 
Figure 11. - Jet centerline shift. 
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Figure 12. - Jet spreading at angle of attack 










